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Abstract 

We compute free energy of quark matter at asymptotically high baryon number density in the 
presence of non-zero strange quark mass including dynamics of pseudo Nambu-Goldstone bosons 
due to chiral symmetry breaking, extending previously existing analysis based on perturbative 
expansion in mg/4/iA. We demonstrate that the CFLiC'' state has lower free energy than the 
symmetric CFL state for < m^/4/iA < 2/3. We also calculate the spectrum of the fermionic 
quasiparticle excitations about the kaon condensed ground state in the regime m^/4/iA ~ 1 and 
find that (m^/4^A)crjt = 2/3 for the CFL-gCFL phase transition, the leading order result reported 
in [ill, is not modified. 
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I. INTRODUCTION 



It is known that while Color-Flavor Locked (CFL) phase is the ground state of 3-color, 

nnn 

3-flavor massless QCD at asymptotically large baryon number density P, y, |^ , the ground 
state is modified by non-zero quark masses. An investigation of this issue has recently 
received a great deal of attention. See reviews and references therein. 

In the CFL phase with zero quark masses, quarks of all three fiavors participate in 
the BCS pairing due to the attractive one gluon exchange potential in the 3 channel. The 
resulting condensate spontaneously breaks the original U (1) b >^ SU (3) coior>^ SU (3) lX SU {3)ji 
symmetry down to the diagonal subgroup, SU{3)coior+L+R, causing the gauge bosons of 
the original SU{3)coior group to become massive by the Higgs mechanism 3^. A linear 
combination of gluon and the photon remains massless; this unbroken gauge symmetry is 
referred to as U{1)q, a "modified" electromagnetism. The low energy bosonic excitations 
are parametrized by S = ^^/^-^^'^'^ /^f^' aj^d B = e^l^lfB (with vr = vr^T") where tt^'s are 
the pseudo-scalar octet of Nambu-Goldstone bosons (NGB's) which arise from the breaking 
of chiral symmetry and T^'s are the fundamental representation SU{3) generators; r/' and (3 
are NGB's due to breaking of the U{1)a (assumed to be a valid symmetry of the theory at 
high density) and U{1)b symmetries, respectively. The CFL symmetry breaking pattern is 
identical to that of QCD at low density j9| with the important difference between the two 
cases being that at high density the system is weakly coupled and the coefficients of the low 
energy EFT may be calculated from the underlying microscopic theory. 

In the presence of non-zero quark masses the ground state changes in the following way. 
A relatively large strange quark mass encourages the system to reduce the strange quark 
number density relative to the density of up- and down-type quarks. The system responds to 
this stress by forming a condensate of collective excitations with quantum numbers of kaons 



in the ground state |lG , 
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12| |. Kaon condensation allows the strange quark number densit 



isity 
0- 



to be decreased without the costly breaking of Cooper pairs in the CFL background 
The ml/ 2 fi term, with ^ being the quark number chemical potential, acts as a dynamically 
generated chemical potential for strangeness; when it exceeds the kaon mass the Bose- 
Einstein condensation of K'^ is triggered and the CFL ground state undergoes a second order 
phase transition to the CFLi^^ phase where the parity and hypercharge are spontaneously 
broken by the condensate. Recently it was argued that taking into account (9(as^s) terms 
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in the effective lagrangian leads to more complicated phase diagram However, in this 
work we neglect this effect for the sake of simplicity. 

So far the analysis has been performed to the leading order in M^//iA, where M = 
diag{mu, rrid, rris) is the quark mass matrix and A is the CFL gap parameter. It was shown 
that in the case of maximal kaon condensate (which occurs if we set m„ = = and 
rris 7^ 0) the free energy was lowered by 5n = -fl (m2/2/i)' /2 oc ml while the hypercharge 
density in the CFhK^ ground state was riy = fl {ml/2 u) , where = (21—8 Log 2)/i^/36 vr^ 
is the leading order result for the pion decay constant Il2l . Il4|. Extrapolating into the 
regime m^/A A ~ 1 one may expect that the kaon condensation free energy is of the same 
order as the BCS condensation energy itself and the number of strange quarks involved 
in the condensate is of the same order as SyU^A/vr^, the (approximate) number of strange 
quarks involved in the BCS pairing. Therefore, it may be important to take into account 
dynamics of NGB's when one studies the properties of the ground state at large m^. 

Recently several calculations of dense quark matter free energy as a function of 



have been performed and C FL-g apless CFL (gCFL) phase transition was observed at 
(m^/4/i A)crit = 1/2 [3, ll^ \v\ . In the gCFL phase the quark-quark pairing still in- 
volves all colors and flavors, but the spectrum contains massless fermionic modes. Similar 
phases were discussed in the context of two flavor QCD ^1^. The objective of this work is 
to obtain better understanding of the CFL ground state as a function of which should 
help better understanding of the CFL-gCFL phase transition and of the properties of gCFL 
phase. 

In Section 2 of this article we calculate free energy of the CFLi^^ phase to all orders 
in the m^/4/iA expansion and show that the CFLif° persists as the energetically favored 
state even for m^/4/iA ~ 1. In Section 3 we extend the leading order calculation of the 
CFL phase quasiparticle excitation spectrum about the CFLK'^ ground state and observe 
that the result reported in Q is not modified by the higher order corrections. 



II. FREE ENERGY CALCULATION 



Let us first list some simplifying assumptions and approximations we make in this work 

• We neglect small light quark masses and assume M = diag{0,0,ms). We work in the 
regime ml/iJ? <^ 1, m^/4A/i < 1 so that the light-strange pairing is still possible 
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jigl - This means that we neglect deviation of Vg = = ^JT^-rn^JJj? , the strange 
quark Fermi velocity, from 1. Note that in this regime the chiral expansion, that is 
an expansion in 9/A ~ m^/AyU, where one retains only few leading terms to attain 
desired accuracy, is expected to break down. In this regime the relevant effective theory 
should include both bosonic excitations of diquark condensate and quasiparticles as 
the degrees of freedom. We will use mean field approximation in which we will freeze 
the "radial" oscillations of the condensate magnitude while allowing oscillations of 
the condensate around the CFL vacuum state corresponding to the NGB's due to the 
chiral symmetry breaking in the massless case. 

Being in the weak coupling regime we will work to the leading order in a^. Also we 
will neglect meson mass terms generated by quark masses. Quark mass terms connect 
particles and antiparticles and the resulting meson mass terms are generally suppressed 
by powers of A/fi oc Exp{— const /qs) and/or M, 2^ - 



• We assume that is small enough so that we are in the CFL phase; we consider only 
color-flavor antisymmetric CFL gap term. 

For excitation energies below fi the relevant degrees of freedom are the nonet of quasipar- 
ticles and holes, ten pseudo NGB's due to the spontaneous symmetry breaking of global 
symmetries and soft gauge bosons. The corresponding effective theory is called the High 



Density Effective Theory (HDET) IjJ, |2l|, |22, |2^. Let's consider the leading terms of the 



HDET Lagrangian of the CFL phase of high density QCD in the mean field approximation. 
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Lj{iv.d - {ml/'iij) - jj,^ - e v.Aem Q)L^ - 9s Tr (^LJL^v.A^^ + 



+ Tr 

+ y e,,keabc X'^'V^.^Ll'p e"'^ + /i.e. + (L ^ i?, X ^ F, /if ^ /if) - 
- h:rG^''G,,--^F^^''F^,. (1) 
Let us define all the quantities in 

• The first term is the mean field potential with coupling constant G being the coupling 
of the effective four fermion interaction generated by the hard gluon exchange. 

• Fermionic fields L\^°'{x) and R\j"'{x) are the left- and right handed quasiparticle fields 
that annihilate states with Fermi velocity = (1,'0), v.v = 1; fiavor indices i,j,k 



take values u,d,s; color indices a,b,c take values r,g,b; spinor indices a,P = {1,2} 
are explicitly shown only in the gap term. Under the original symmetry group SU (3)c x 
SU{3)l X SU{3)r, transforms as (3,3, 1) and transforms as (3, 1,3). 

We separated the SU(3) singlet part of the leading HDET mass term MV2/i. The 

singlet term will only be responsible for the baryon density shift due to non-zero 

quark masses, but will not affect the dependence of the condensation energy and 

from now on will be neglected. The remaining mass dependent terms are 

, MMt 1 MMt 

= Tr , 

2fi 3 2/i ' 

„ MtM 1 MtM 

uf = Tr . (2) 

'^^ 2/i 3 2/i ^ ^ 

With the M we use fx^ = fxf = —fXs Tg, where /t^ = ml/\^fi. 

The NGB fields are collected in X'^^ ~ e^'^ e'^'^Ll^Llj) and V"" ~ e'''^ e""'' {RliRlj) , 
3x3 unitary matrix valued composite fields describing oscillations of the left and 
right handed quark-quark condensates about the CFL ground state. Under the original 
symmetry group5'[/(3)c x SU (3) l x SU{3) r, X transforms as (3, 3, 1) and Y transforms 



as (3,1,3) |lj,|24|. 



• G^'^ and Ff^^ are the usual field strength tensors for the gauge fields and ^em- 

So, in our approximation the strange quark mass only appears in the role of an effective 
chemical potential for strangeness. An infinitesimal induces kaon condensation in the 
ground state. What changes as increases? As was stated above, the kaon mass terms 
generated by the quark masses are always suppressed by the powers of A//i and/or as 
and cannot counterbalance the "chemical potential", m1/2fi, driving the Bose- Einstein 
condensation. For as long as the CFL pairing pattern is preferred, there is nothing that 
can significantly alter the condensate in the ground state as increases. And in our 
approximation there is literally nothing that can change the form of the kaon condensate in 
the ground state. This is what our calculation will demonstrate. 

Let us start by changing basis of quasiparticle fields to diagonalize the gap term. We 
write 



A=l 



(3) 



where for A = 1, . . . , 8 are the Gell Mann matrices and Ag" = y2/35*" The gap 
term becomes 



Cgap = ^Bet[X]AABL^,aL^/3 + h.c. -{L^R,X^Y) 



(4) 



where Aab = ^AiX^/X^^ - A^^A^) = -A,A,B = 1..8 and for A,B = 9, Aab = |A(| x 
3 X 3 - I X 3) = 2A and we used e^tc^'^^X^^X^'^ = Det[X]e™"^ In this calculation we 
neglect Ub{1) and Ua{1) NGB's and set Det[X] = Det[Y] = 1. The Lagrangian ^ may be 
rewritten as 



£ = -2 X + 



G 



At 



'AB 



A,B=1 



At 



-A,5^B 



\ 



-AaS 



AB 



where v = (1, —-0) and p = v.p — The rest of the definitions is as follows: 



\ 



and 



R 



( 



{R 



V 



are the components of the Nambu-Gor'kov field; 



(5) 



(6) 



(7) 



AB 



-Tr\AX\iv.d - /if - ev.AemQ)XXB - -gsTr XaXbv.A 



T 
c 1 



(8) 



AB 



-Tr XAY^iv.d - /if - e v.A^^ Q)YXb - - (?s Tr A^As v.A] 



(9) 



Before we proceed let us state the additional simplifying assumptions we make. 
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We will neglect oscillations of the NGB fields about their ground state expectation 
values as the corresponding contribution to the free energy density is suppressed by 
(A//i)^ relative to the fermionic contribution. The values of X and Y in the ground 
state are to be determined by solving the corresponding equations of motion. 

We will also neglect oscillations of the gauge fields about their ground state expecta- 
tion values. The (A)'s are to be determined by solving the corresponding equations 
of motion which in the case of homogeneous quark matter are also the conditions for 
the gauge charge neutrality. These oscillations are important for the calculation of 
the gap parameter as the gap equation is dominated by the soft magnetic gluons j^l j 
but here we do not attempt the calculation of the gap and neglect the corresponding 
terms in the effective action as being of higher order in the coupling expansion. Our 
approximation is then equivalent to working in the model with a color-flavor anti^m- 
metric short range interactions in the mean field approximation (N JL model) [2, ■ 
In our approximation the only role of the gauge fields is to ensure the gauge charge 



neutrality of the ground state 
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27l |. The spatial components of gauge fields do not 
have sources and, therefore, vanish in the ground state. Then the kinetic terms for 
the gauge fields vanish, too. 

So, in this approximation we consider dynamics of quasiparticles in some gauge charge 
neutral meson field background. To calculate free energy we need to determine what this 
background is. 

Starting from (0) let us formally integrate out fermions and then expand in powers of 
X and y. This is the usual EFT expansion; effectively we expand in m^/4/iA< 1. Let us 
stress that we do not truncate this series at any finite order, that is, we do not rely on the 
smallness of the expansion parameter. We do rely on the analyticity of the free energy as a 
function of m^; our analysis is valid within the radius of convergence of the expansion and 
breaks down at the point of a phase transition. We get the following effective action 

S[X,Y,A] 
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+ - y — Tr 



S 



AB 

V 



s 



AB 

V 



BA 



-3^^; 




2n 



(10) 



where 



X.^^ = -Tr \a [X\iv.d - - t;. A„ Q)X - v.A^ 



(11) 



and 



= ^Tr [Ft(,^;.5 - /if - t;.Ae^g)r - v.A 



(12) 



and where S ^ is the inverse massless CFL propagator 



(13) 



Also we have rescaled gauge fields to eliminate explicit dependence on the couplings. The 
operation Tr in |TO|l means trace over the 18 dimensional color fiavor Nambu-Gor'kov space 
as well as over the space-time indices. The leading order terms in the expansion (fTTH) produce 



leading order terms of the 



spatially uniform fields) 
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ow energy EFT Lagrangian (for simplicity here we consider 



28| 



-YriX^D^Xf + llr{j\D^Yf\ + O (^)M , 



(14) 



where the covariant derivatives are given by 



D^X = doX + iXAf + iAl'^QX + i/if X, 
D^Y = doY + lYAf + lA^'^QY + «/if y. 



(15) 



For spatially homogeneous time independent configurations neutral under \Jq{\\ i.e. for 
[X, Q] = and [F, Q\ = 0, it's easy to see that 



A = Af + A^J^Q = -i(Xt/ifX + Y^fxfY), 



(16) 



the leading order solution, satisfies the equations of motion 6S/6 Aq = for any term in the 
expansion (fTUI) [27]. Since the CFLK^ is devoid of electrons, we set A^"^ = 0. If we constrain 
the ground state to be the symmetric CFL configuration and set X = 1, y = 1, then 



A = Af + Al"" Q 



(17) 



and Aem = 0, = 0, 

The fact that the leading order solution A satisfies equations of motion at any order 
in the expansion follows from the general form of gauge invariant effective action of a 
superconducting system [2^. 

Let us integrate out color gauge fields which in this approximation is achieved by simply 
substituting (fT^ back into the effective action. Now we have 



Sin 



3|A| 



Tr Log fS"^ + Jy + Ja] + Jv - -l 



(18) 



where 



Jv 



TrA^[(Xt/ifX + FVfF),AB] 



J 



V 



-TrAB[(XVf'X + FVfF),A^] ) 







A 



J 



(19) 



(20) 



4 1 -Tr AB(Xt/ifX - Y^fi^Y)XA 

Considering only spatially homogeneous time independent configurations, we expand ()18|) 



in m1/4fiA and obtain 
= - 



o|A|2 00 1 

2 I -L^-zTrLogS-'+tJ2 ^Tr [S S J^] 

^ n=l ^™ 



2n 



(21) 



Here Tr also means trace over the 18 dimensional color fiavor Nambu-Gor'kov space as well 
as over the space-time indices. The leading term in the expansion corresponds to the leading 
term of the CFL xPT Lagrangian 



J n 



.2 A 2 / "^s \4 



(22) 



where S = XY'' and Dq S = ifi^ S — fif, the Bedaque-Schaefer covariant derivative for 



a constant field 
may show that 
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. We observe that Jy terms are not present in the expansion. One 



E — T,KO — 



1 o^ 

i 
i J 



(23) 



the leading order solution for the meson field |ll|, Q , satisfies equations of motion for E 
for any term in ()21|) . From S we may determine values of matrices X and Y in the ground 




FIG. 1: Free energy difference 6Q = = S^o) ~ = 1))/I^c'FlI' between S = S^^o and 

S = 1 ground states as a function of m'^/AiJ.A. The massless CFL BCS condensation energy, Q^pj^, 
is defined in The ii'*^ condensed state has lower value of free energy. 



state which in the unitary gauge are given by X^o = ^i^o, Yko = [2^ 

Now that we know the background gauge field values that ensure gauge charge neutrality 
and fields X and Y in the ground state, we may calculate the difference in free energies 
between CFLK^ and the SU (3) symmetric CFL ground state. The expression is 



6Q 



3|AP 



AS= 11 



2 X 



2 X 



r 9 

— — + 2 vr / ^ ^s) - >^a{p; m^, a = a, S = S^^o 

3|AP /■ ^ 

2 / E (A!i(p; m,) - Aa(p; m„ A = i, S = 1 



G 



'pa=i 



(24) 



where {pO — Aa(p; m^, A, S), pO + \A{p',Tns, A,I])} are the eigenvalues of the full inverse 
propagator in and {pO — A° (p;ms), pO + \^{p;ms)} are the eigenvalues for the charge 
neutral normal phase. Note that we don't need an explicit expression for the latter. We 
neglect the difference in the value of the gap parameter for the two cases which we have 
checked to be small. Alternatively, we calculated the difference in free energy starting from 
fll8|l . The results are shown in Fig. d The condensed configuration has lower free energy 
relative to the symmetric CFL state even for m^/4/iA ~ 1. The result of hypercharge density 
calculation in the CFLii'^ phase is shown in Fig. |21 We note that the numerical results are 
surprisingly close to the leading order EFT predictions. 
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ny/nBCS 




0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 2: The hypercharge density ny/nBCS in the CFLK^ ground state as a function of m^/4/iA, 
where ubcs = 3^^ A/vr^, the approximate concentration of quarks of a given flavor involved in 
the BCS pairing. The result is close to the leading order EFT prediction. 



In this calculation we performed the free energy renormalization by trading the bare 
coupling 1/G for the (physical) superconducting gap parameter Aq of the massless CFL by 
expressing 1/G from the gap equation 



sn 

6A 



0. 



(25) 



ms=0,A=Ao 

Then, for example, the massless CFL BCS condensation free energy is given by 



^CFL = -3 /i' AI/tt' 



(26) 



III. QUASIPARTICLE SPECTRUM 

Based on the results obtained, in this Section we extend the leading order calculation of 
the ,uas,pa.io.e spect.u. ,n the p.esenee of ka„„ e„nde„sed bae.gtou„d perfo^ed ,n Q 

into the regime m^/4/iA ~ 1. Starting from (fTHj) we write down the leading terms of the 
effective Lagrangian for the nonet of quasiparticle excitations about CFLA'° ground state 

C = Tr'iljljv.D'iljLv + Tr'i/jl^Ao'ipLv + ^^i't-vi'Lv + h-c. + 

TrijjliJv.D^R^ - TrT/'Jj^A - ^^^i^v^Rv + h-c, (27) 
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where the nonet field has been defined in Q- This Lagrangian may be rewritten in the 
form prescribed by the chiral perturbation theory for the baryon fields js^ 

C = Trtpliv.Dtp^ + DTriljl'y5{v.A,ipv} + FTriljl'y5[v.A,ip^] + 

+ %^nC'75^f + /^.c., (28) 

with D = F = 1/2 and C being the charge conjugation operator. The covariant derivative 
here is D^^ ip = 9^ V^ + ifV^, the time components of the vector and axial currents modified 
by the quark mass terms are 

Vo = ^ {do + ^^i^)C + a^o + ^ /if) = (^Vf ^ + Cf^f^) , (29) 

and 

A = ^ (do + {do + ^^if) e) = -\ (^Vs^ - ^/if^^) • (30) 

Shown in Fig. El are the masses of the quasiparticles that become gapless as a function of 
m^/AfiA. The S = 1 case develops gapless modes at rn^/AfiA = 1/2, while a CFLi^''^ mode 
becomes gapless at m^/4/iA = 2/3. A detailed description is presented in Here we only 
observe that the leading order in m^/4/xA result and the result of the full calculation coincide. 
It follows from the fact that in the approximation we use the leading order in ml/ An A 
expressions for the background gauge fields and S are not modified by the higher order terms. 
It has been argued that the appearance of a gapless mode signifies the point of CFL-gCFL 
phase transition ^|. Our calculation confirms that due to the nontrivial fiavor alignment 
in the ground state the CFL-gCFL phase transition is expected at (m^/4yuA)c = 2/3 which 
is 4/3 times higher than the symmetric CFL value. 

Note that a non analyticity at m^''** is built into the expression for free energy ()24j] since at 
rris = rrig"^ the momentum integral over > has to be modified to take into account the 



breach in the quasiparticle spectrum. Then the value of mf^^ = y8/iA/3is the convergence 
radius of the series expansions (fTUI) . (PT|) which is consistent with the fact that at this 
point the system undergoes a phase transition. 

We have neglected the quantum mechanical oscillations of the collective bosonic fields 
about CFLK^ ground state. One may have to worry about the stability of the ground 
state since as grows several quasiparticles are becoming lighter, whereas the massive 
mesonic excitations about CFLi^'" ground state have masses of order m^/2/x and at some 
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FIG. 3: Masses of the quasiparticles that are lowered due to the presence of function 
of m'l/AfiA. The dashed hne depicts a doubly degenerate E = 1 case mode, solid lines are the 
S = T,j^o case modes. The results are in agreement with Q]. 

large enough rris a heavy Goldstone may decay into a particle-hole pair. However, due to 
the hypercharge and isospin charge conservation the quasiparticle-hole pair in the fermionic 
loop may not go on-shell for any rris and the ground state is stable. 



IV. DISCUSSION: SOME OPEN QUESTIONS 

We have demonstrated that taking into account dynamics of pseudo NGB's in the CFL 
phase of high density QCD and allowing a nontrivial flavor alignment in the ground state 
gives a sizable decrease in the free energy in the regime / A ~ 1 . It is remarkable that we 
were able to find a sensible approximation where it was possible to do the calculation to all 
orders in the rris expansion. We also extended the leading order calculation of quasiparticle 
spectrum modified by the presence of kaon condensed background into the regime rn^/AfxA ~ 
1 and found the leading order result unmodified. 

Several comments may be made based on the results of this work. First, since the CFL- 
gCFL phase transition is continuous, the kaon condensate in the gCFL ground state cannot 
be turned off abruptly and will persist at least for some values of rrig > m^"*. While we argued 
that in the CFL nothing can significantly alter the magnitude of the kaon condensate, the 
non-zero non CFL gap terms, Ai and A3 defined in ()31|) . may just do that by generating 
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mass terms for the kaons. The question whether condensate will persist in the gCFL 
phase requires further investigation. 

The most interesting issue that we have not addressed in this work is the investigation of 
the CFL-gCFL phase transition and of the gCFL ground state including dynamics of pseudo 
NGB's. In such a calculation one should allow deviations from the CFL pairing pattern in 
the gap term. For example, one may consider 

A"* = AcflS"" + A^^3^ (31) 

with 

Ar3'= = d^a(7(Al,0,A3). (32) 

This pairing ansatz introduced by Alford et al in Q| is the simplest one that allows both 
disruption of the light-strange pairing (A3 7^ 0) which is directly induced by m^, and of the 
pairing between light quarks (Ai 7^ 0) which is due to the gauge charge neutrality and which, 
of course, is also indirectly induced by the strange quark mass. Unlike the CFLif*^ case, in 
the presence of the non CFL gap terms one does not know anymore the exact values of the 
background gauge fields and S in the ground state. A perturbation theory in Ai 3 about 
the CYhK^ state would seem to be a reasonable approach. To perform this calculation one 
needs to construct a low energy EFT in the regime m^/4/iA ~ 1 describing dynamics of 
light bosonic and fermionic modes. 

In this simplified calculation we treated gluons as classical background fields. It would be 
interesting to repeat the calculation of the CFL free energy including the effects of quantum 
mechanical oscillations of the gluon fields into account. 

These questions are left for future work. The full answer to the question: "What is the 
ground state of dense QCD as a function of quark masses?" is yet to come. 
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